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ABSTRACT 


A  computer  simulation  has  been  used  to  calculate  the  sampling 
efficiency  of  a  downward-facing  aerosol  sampler  having  a  large  entrance 
aperture.  Flow-field  calculations  have  been  done  by  solving  the 
Laplace  equation  for  a  flat  disk  with  uniform  air  velocity  over  its 
surface.  Initial  positions  from  which  particles  reach  an  equilibrium 
position  at  which  the  sedimentation  speed  equals  the  air  suction  speed 
of  the  sampler  are  obtained  by  iterative  solution  of  the  equations  of 
motion.  These  distances  are  time-independent  and  are  called  critical 
distances.  The  time-dependent  sampling  efficiency  is  obtained  by  com¬ 
bining  critical  distances  with  the  sampler  interaction  distance  for  a 
given  sample  volume  of  air.  The  effect  on  sampling  efficiency  of 
various  parameters  such  as  the  sampling  flow  rate,  sampling  period, 
height  of  the  aerosol  above  the  filter,  and  rim  size  is  evaluated. 

RESUME 


Une  simulation  sur  ordinateur  a  6te  utilisee  pour  calculer  1'ef¬ 
ficacite  d'un  echant i 1 lonneur  d'aerosol  inverse  S  grande  ouverture. 

Les  champs  d ' ecouleraent  ont  ete  calcules  en  solutionnant  l'equation  de 
Laplace  pour  un  disque  plat  dont  la  vitesse  de  l'air  est  uniforme  sur 
la  surface.  Les  positions  initiales  desquelles  les  particules  aboutis- 
sent  en  equilibre  ^  une  position  ou  la  vitesse  de  sedimentation  egale 
la  vitesse  de  succion  de  l'air  de  1 ' £chanti llonneur  sont  etablies  par 
Iteration  des  equations  du  mouveraent.  Ces  distances,  appelees  dis¬ 
tances  critiques,  sont  independantes  du  temps.  La  dependance  tempo- 
relle  de  1'efficacite  d' echantillonnage  est  obtenue  en  combinant  les 
distances  critiques  et  la  distance  d ' interaction  pour  un  certain  volume 
d'air.  L'influence  sur  1'efficacite  d ' echant 1 1 lonnage  des  parametres 
tels  que  le  debit  et  la  periode  d ' echant i llonnage,  la  hauteur  d'aerosol 
au-dessus  de  1' echant i 1 lonneur  et  la  largeur  du  rebord  est  calculee. 


j'Kibn.jy  Cedes 

n.r*  i  Avj'(  J  ri/c>r~ 
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NOMENCLATURF 


a 

A 

AC 

ST 

AS 

B 

BC 

F, 

8 

h 

H  (s,T) 

Jl 

L 

Q 

r 

R 

Ri 

s 

S  (s,T) 


sampler  aperture 

semi-minor  axis  along  R  (representation  of  a  potential  line 
by  an  ellipse) 

semi-major  axis  along  R  (representation  of  critical  position 
AS  by  an  ellipse) 

mean  interactive  distance  during  the  sampling  period  T 

critical  position  along  R  for  a  given  and  s 

serai-major  axis  along  Z  (representation  of  a  potential  line 
by  an  ellipse) 

semi-minor  axis  along  Z  (representation  of  critical  position 
AS  by  an  ellipse) 

sampling  efficiency 

gravitation  constant  981  cm/s2 

height  of  aerosol  above  the  sampler 

function  related  to  the  height  of  the  aerosol  above  the 
sampler 

Bessel  function  of  order  1 
sampler  radius 
sampling  flow  rate 
particle  radius 

particle  coordinate  along  R-axis 
initial  particle  coordinate  along  R-axis 
pr2  -  defines  particle  mass  size 

function  related  to  the  surface  crossed  by  particles  under 
sedimentation  and  to  the  sampled  volume  of  air 


T 


sampling  period 
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time  required  for  the  air  at  position  ACCsj)  to  reach  the 
filter  paper  surface 

time  required  for  the  air  at  position  L  to  reach  the  filter 
paper  surface 

volume  of  air  sampled  during  period  T 

sample  volume  of  particles  of  type  s  from  above  the  filter 
sample  volume  of  particles  of  type  s  from  below  the  filter 
volume  of  air  sampled  during  period  ^ 

volume  of  air  sampled  during  period 

R-component  of  air  velocity  in  air  flow  created  by  the 
filter  holder 

settling  velocity  of  a  particle  of  parameter  s 

Z-component  of  air  velocity  in  air  flow 

particle  coordinate  along  Z-axis 

position  along  Z-axis  where  VZ(Z,  0)  =  j V g J 

initial  particle  coordinate  along  Z-axis 

initial  coordinate  value  along  Z-axis 

air  viscosity 

current  line  function 

particle  density 
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1.0  INTRODUCTION 


Aerosol  sampling  has  been  the  object  of  much  investigation.  The 
purpose  of  aerosol  sampling  is  to  determine  the  concentration  and  the 
mass  or  size  distribution  of  particles  suspended  in  a  given  volume  of 
air  at  a  given  time.  A  fundamental  parameter  required  by  the  re¬ 
searcher  is  the  sampling  efficiency,  defined  as  the  fraction  of  parti¬ 
cles  of  a  given  size  and  density  that  will  be  collected  by  a  sampling 
system. 


The  sampling  efficiency  changes  from  one  sampling  system  config¬ 
uration  to  another.  Calculations  usually  deal  with  the  sampler  orien¬ 
tation  and  air  flow  rate.  Non-lsokinetic  conditions  as  studied  by  W.M. 
ter  Kuile  (Ref.  1)  are  discussed  in  his  paper  entitled  "Comparable  Dust 
Sampling  at  the  Workplace  (a  contribution  to  the  standardization  of 
"total  dust"  measurements)".  He  draws  the  conclusion  that  the  downward 
inlet  arrangement  best  meets  the  requirements  for  standardization.  His 
conclusion  was  based  on  the  work  of  Walton  (Ref.  2),  Levin  (Ref.  3), 
and  Davies  (Ref.  4)  who  developed  theoretical  models  based  on  a  point- 
source  sink  with  a  thin-walled  filter  holder.  Parameters  such  as  the 
sampling  period  and  the  height  of  the  aerosol  above  the  filter  were  not 
considered. 

The  present  work  involves  a  computer  simulation  of  the  inter¬ 
action  of  particles  undergoing  sedimentation  on  a  large-aperture  in¬ 
verted  filter  holder.  Flow-field  calculations  are  done  by  solving  the 
Laplace  equation  for  a  flat  disk  with  uniform  air  velocity  over  its 
surface.  Critical  distances  are  then  obtained  by  iterative  solution  of 
the  equations  of  motion.  These  distances  are  time-independent.  The 
time-dependent  sampling  efficiency  is  obtained  by  combining  critical 
distances  and  the  sampler  interaction  distance  for  a  given  sample  vol¬ 
ume.  Thus,  the  effect  of  various  parameters  such  as  the  sampling  flow 
rate,  sampling  period,  height  of  the  aerosol  above  the  filter,  and  rim 
size  can  be  evaluated.  The  sampling  setup  is  illustrated  in  Fig.  1. 
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This  work  was  conducted  at  DREV  between  February  and 
December  1983  under  PCN  2IB05,  Aerosols. 

2.0  DETERMINATION  OF  CRITICAL  POSITIONS 

The  critical  position  of  a  particle  of  radius  r  and  density 
is  the  closest  position  to  the  sampler  at  which  the  particle  will  not 
be  collected  by  the  sampler  and  thus  will  escape  it.  As  long  as  the 
settling  velocity  of  a  particle  (Vg)  is  lower  than  the  airspeed  on  the 
filter  paper  of  an  inverted  sampler,  critical  positions  will  exist. 

In  order  to  calculate  critical  positions,  the  equations  of  mo¬ 
tion  of  a  particle  in  a  flow  field  created  by  an  inverted  filter  holder 
have  to  be  solved. 

Figure  2  illustrates  the  interaction  of  particles  with  the  flow 
field  created  by  the  inverted  filter  holder,  the  radius  of  which  is 
a  =  6.7  cm.  Particle  trajectories  as  indicated  are  only  examples  but 
the  current  lines  shown  are  those  created  by  a  flat  disk  with  uniform 
air  velocity  at  Its  surface.  It  has  been  assumed  that  settling  parti¬ 
cles  start  Interacting  with  the  field  only  for  Z  <  0.  The  middle  tra¬ 
jectory  defines  two  critical  positions:  at  Z  »  0,  R  «=  AC  and  at 
Z  -  BC,  R  -  0. 
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P 


(vr(z,r)  -  -£) 


[2] 


where 


g:  981  cm/s2 

r:  particle  radius  (cm) 

t :  t ime  ( s ) 

R:  particle  coordinate  along  R-axis 

Z:  particle  coordinate  along  Z-axis 

V Z-component  of  air  velocity  in  air  flow  created  by  the  filter 
holder 

V  :  R-component  of  air  velocity  in  air  flow  created  by  the  filter 

K 

holder 

y :  air  viscosity 

p  :  density  of  the  particle. 

P 

Velocities  V  and  V  are  calculated  using  the  following 
iu  K 

expressions: 


VZ(Z,R) 


-l  H(Z,R) 

r  snr 


[3] 


vr(z,R) 


1  3<KZ,R) 

\ - TT~ 


[4] 


where 


<KZ,  R) 


zl 

it  a 


00  kz 

/  J^kR)  Jj  (ka)  dk 

o 


[5] 


where  Z  <  0 
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Q  :  sampling  flow  rate  (SFR)  (cm3/s) 

a  :  radius  of  the  circular  aperture 

Jj(  ):  Bessel  function  of  order  1 
k  :  variable  in  Fourier-Bessel  space. 

Refer  to  Appendix  A  for  relevant  details. 

Prior  to  the  solution  of  the  equations  of  motion,  it  is  essen¬ 
tial  that  the  velocities  V„  and  V  be  calculated.  In  order  to  do  so, 
the  function  is  calculated  over  the  area  of  interest  for  m  x  n  posi¬ 
tions.  Using  the  definition  of  the  derivative  for  the  calculation  of 

velocities  V  and  V  ,  two  matrices  of  size  (m-1)  x  (n-1)  are  defined. 

Z  JK 

Then,  when  the  numerical  solution  of  the  equations  of  motion  are  re¬ 
quired,  values  of  flow  velocity  are  obtained  by  interpolation  (see 
Appendix  B  for  more  details). 

2.2  Numerical  Solution  of  the  Equations  of  Motion,  and  the  Criterion 
for  Collection _ _____ _ 

The  equations  of  motion  are  numerically  solved  using  the  fourth- 
order  Runge-Kutta  technique.  Initial  conditions  of  the  position  and 
velocities  are  as  follows: 


[6] 
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The  equations  of  motion  are  solved  for  a  given  particle  size  and 
density.  However,  we  will  refer  to  a  more  general  quantity,  s,  defined 
as  p^r2.  Particles  of  equal  s  will  act  in  the  same  way. 

The  criterion  for  collection  is  simple  because  of  the  azimuthal 

symmetry  of  the  problem:  a  particle  will  be  collected  if,  at  any  time, 
d  z 

-pC-  becomes  greater  than  zero.  A  particle  will  not  be  collected  if  at 

any  time  Z  <  (Z  is  negative-valued)  where  is  the  position  at 

which  V_ (Z  ,  0)  =  V  .  This  is  on  the  axis  where  the  velocity  is 

Z  c  Is! 

maximum  in  the  Z-direction. 


Once  the  status  of  a  given  particle  is  known  (whether  it  will  be 


collected  or  not),  the  initial  value  is  increased  or  decreased.  Its 
status  is  redefined  and  so  on  until  the  difference  between  the  initial 
position  of  a  collected  particle  and  of  an  uncollected  particle  is 
smaller  than  a  predetermined  accuracy  AR  ( -  0.1  cm).  This  position  is 
called  the  critical  distance  AS  for  a  given  Z^  and  s.  A  particle  with 
these  initial  conditions  Is  brought  to  R  =  0  and  stays  motionless  in 


the  air  at  a  distance  Z  where  VZ(Z,  0) 


The  particle  weight  and 


the  air  stream  drag  force  are  then  in  equilibrium. 


2 .3  Representation  of  Critical  Positions  on  the  Z-R  Plane 

Critical  positions  were  calculated  for  a  range  of  values  of  the 
parameter  s  at  different  initial  heights  (Z^)  for  the  two  sampling  flow 
rates  Q  of  30  L/min  and  50  L/min.  Figure  3  shows  critical  positions 
for  an  s  of  90  and  160  g  um2/cm3;  the  continuous  curves  are  potential 
lines  created  by  a  disk  of  radius  6.7  cm  with  uniform  velocity  over  its 
surface.  Critical  positions  can  be  represented  by  a  family  of  ellipses 
with  principal  axes  AC  and  BC  along  the  R  and  Z  axes  respectively. 
Figures  4  and  5  show  critical  positions  of  the  principal  axes  AC  (+) 
and  BC  (-)  as  a  function  of  s  for  sampling  flow  rates  of  30  and 
50  L/min.  The  smaller  the  particles,  the  farther  the  critical 
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position.  AC  is  always  greater  than  BC.  Figure  6  shows  critical  posi¬ 
tion  AC(s)  as  a  function  of  BC(s).  Sampling  flow  rate  data  for 
30  L/min  are  represented  by  (•+;  while  those  obtained  at  50  L/min  are 
represented  by  (-)•  The  two  curves  obtained  agree  fairly  well  with 
each  other.  For  large  values  of  BC(s),  it  appears  that  there  is  a 
linear  relationship  between  the  two  principal  critical  position  axes: 

AC(s)  =  1.6  BC(s)  [ 7] 

This  expression  may  be  used  for  extrapolation  of  critical  posi¬ 
tions  for  small  s.  Since  BC(s)  is  easily  calculated  by  finding  the 
position  on  the  Z-axis  where  VZ(Z,  0)  «  Vg,  critical  positions  can  be 
obtained  easily  using  [7]. 

3.0  TIME-DEPENDENCE  OF  SAMPLED  VOLUME  OF  AEROSOL 


Critical  positions  have  been  established  for  time-independent 
sampling.  The  volume  of  air  QT  (Q:  sampling  flow  rate,  T:  sampling 
period)  drawn  through  the  filter  defines  the  space  where  the  air  has 
been  taken  from.  This  space  is  bounded  by  a  potential  line  which 
crosses  the  Z  and  R  axes  at  positions  B  and  A.  As  for  critical  posi¬ 
tions,  a  good  representation  of  the  potential  line  could  be  given  by  a 
family  of  ellipses.  A  and  B  become  halves  of  the  principal  axes  of  an 
ellipse.  Figure  7  shows  A  and  B  versus  sampled  volume.  The  volume  has 
been  obtained  by  the  evaluation  of  the  integral  for  the  rotation  of  the 
potential  line  around  the  Z-axis.  So,  for  example,  for  a  sampling  flow 
rate  of  Q  ■  30  L/min  and  a  sampling  period  of  T  -  1  min,  a  volume  of 
3.0  x  104  cm3  of  air  is  sampled.  The  rotated  potential  line  (surface) 
that  bounds  this  volume  is  half  an  ellipsoid  of  revolution  with  half¬ 
axis  A  and  B  equal  to  24  and  25  cm. 

Particles  collected  by  the  inverted  sampler  come  from  two  areas: 
particles  below  the  plane  Z  *  0  at  t  -  0,  and  those  above  the  plane 
Z  ■  0  at  t  ■  0.  For  particles  below  the  plane  Z  »  0,  the  volume 
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FIGURE  7 


FI' 


swrx-  V3.VE  (c=!) 

Elliptical  representation  of  sampled  volume.  The  volume 
delimited  by  potential  lines  created  by  the  sampler  are 
represented  by  ellipses  of  revolution  around  Z-axis 


;"RU  8  -  Sampled  volume  of  air  and  critical  positions 
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bounded  by  critical  positions  and  its  interaction  with  the  sampled 
volume  are  considered.  For  particles  above  the  plane  Z  *=  0,  the  height 
of  the  aerosol  above  the  plane,  the  surface  delimited  by  the  critical 
distance  at  Z  *=  0,  i.e.  AC,  and  the  interaction  distance  A  of  the  sam¬ 
pler  are  considered. 


Sampled  Volume  of  Aerosol,  Particle  Below  the  Filter  at  t 


A  sampled  volume  QT  defines  a  volume  delimited  by  a  rotated 
potential  line.  This  potential  line  and  critical  positions  can  be 
represented  by  ellipses  with  axes  A,  B,  AC(s)  and  BC(s). 


Three  different  cases  have  to  be  considered  (see  Fig.  8). 


Case  #1:  The  volume  delimited  by  critical  positions  of  a  parti¬ 
cle  of  given  s  is  completely  inside  the  sampled  volume 
of  air.  The  sampled  volume  of  particles  of  type  s 
V^(s,t)  equals  the  volume  delimited  by  the  critical 
positions. 


Case  112:  The  volume  delimited  by  critical  positions  for  a  given 
s  is  completely  outside  the  sampled  volume  of  air. 

The  sampled  volume  V,(s,t)  of  particles  of  type  s 

D 

equals  the  sampled  volume  of  air. 


Case  H3:  The  volume  delimited  by  critical  positions  for  a  given 

s  is  partially  outside  the  sampled  volume  of  air.  The 

sampled  volume  of  particles  of  type  s  V_(s,t)  equals 

b 

the  volume  of  the  shaded  area. 


The  procedure  used  to  calculate  the  sampled  volume  of  an  aerosol 


of  type  s  is  as  follows: 
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a)  polynomial  representation  of  sampled  volume  of  air;  princi¬ 
pal  axes  A  and  B  as  functions  of  QT; 

b)  polynominal  representation  of  critical  positions;  principal 
axes  AC  and  BC  as  functions  of  s; 

c)  identification  as  case  1,  2  or  3; 

d)  evaluation  of  sampled  volume. 

3 .2  Sampled  Volume  of  Aerosol,  Particle  Above  the  Filter  at  t  ■  0 

A  particle  undergoing  sedimentation  will  interact  with  the  sam¬ 
pled  air  when  It  crosses  the  plane  Z  *  0.  The  sampled  volume  of  parti¬ 
cles  V^Cs.T)  Is  equal  to  H(s,T)  x  S(s,T).  These  two  functions  repre¬ 
sent  the  height  of  the  aerosol  and  the  equivalent  surface  crossed  by  a 
particle  of  type  s  for  a  sampling  period  T.  Figure  9  shows  the  various 
parameters  that  have  to  be  considered  for  the  evaluation  of  these  func¬ 
tions.  The  parameter  h  is  the  height  of  the  aerosol  above  the  sampler, 
the  parameter  a  is  the  radius  of  the  aperture  (a  «  6.7  cm),  L  is  the 
radius  of  the  sampler,  A  is  the  major  axis  of  the  ellipse  delimiting 
the  sampled  volume,  of  air,  and  AC(Sj  )  and  AC(s2 )  are  the  axes  delimi¬ 
ting  these  critical  distances  for  particles  of  type  Sj  and  s2  . 

3.2.1  Evaluation  of  S(s,T) 

The  function  S(s,T)  is  defined  by: 

S(  8  ,T)  -  x  (TT(s,T)  -  L2  )  for  TC  >  L 

*  0  for  TC  <  L  [  8] 

where  TU(b,T)  Is  the  mean  Interactive  distance  during  the  sampling 
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period  T.  Particles  of  type  Sj  and  S2  represent  the  two  cases 
possible. 

Case  #1,  AC(s)  >  A 


Using  the  mean  value  theorem: 


jfC(s,T) 


T 

/  A(t)dt 


[9] 


where  A(t)  Is  a  function  which  represents  the  diminution  of  A  as  sam¬ 
pling  proceeds,  is  the  time  required  for  the  air  at  position  L  to 
reach  the  filter  paper  surface. 


Case  92,  AC(s)  <  A 


^ACf s  }  1  AC(s2 ) 

2T(s,T)  -  (1  -  -pS2;)  AC(s2)  +  (_ - i— )  /  A(0  dt 


AC(s2) 


L  T, 


[10] 


where  T  is  the  time  taken  by  the  air  at  position  AC(s?)  to  read 

AL  I  S2  )  *■ 

the  filter  paper  surface. 


Using  the  expressions  V  *  QT, 
the  equations  for  cases  1  and  2  become: 


QTL  and  VAC(s2) 


QT 


AC(s2)' 


AC(s )  >  A 


AC(s,t)  -  y  l  y  f  A(V)dV 


[11] 
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AC(s )  <  A 


ac(s>T) 


(i 


V 

AC(s2  ) 

- 7 - ) 


AC(s2) 

AC(s2)  +  (_ - i - _)  /  A( V)dV 

AC(e2)  L  V 


[12] 


Figure  6  shows  the  relation  between  A  and  V.  A  polynomial 
representation  of  A  versus  V  is  produced  and  then  the  integration  can 
be  easily  performed. 


3.2.2  Evaluation  of  H(s,T) 

The  function  H(s,T)  represents  the  height  of  the  aerosol  which 
crosses  the  interactive  surface  as  a  function  of  the  parameter  s,  the 
height  h  of  the  aerosol  above  the  filter,  and  the  sampling  period  T. 


The  expression  for  the  settling  velocity  and  the  speed  given  by 
the  ratio  h/T  defines  a  particle  of  parameter  sc  such  that  all  parti¬ 
cles  with  s  greater  than  sc  cross  the  interactive  surface  in  a  period 
of  time  shorter  than  the  sampling  period  T.  For  particles  with  s 
smaller  than  s^,  only  a  fraction  at  height  h  cross  the  interactive 
surface.  This  fraction  is  equal  to  the  ratio  of  settling  velocities  or 
s/s^.  However,  in  both  cases  the  height  of  the  aerosol  crossing  the 
interactive  surface  has  to  be  multiplied  by  the  fraction  of  the  sampled 
VL 

volume:  (1  -  ^— )  . 


The  function  H(s,T)  is: 

V. 

H(s,T)  -  (1  -  yZ.)  •  h  for  s  >  sc  [  13] 

VL 

'  (1  -  ~>  *  h  *  - 


[U] 
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where 


sc  -  9yh/2gT  [  15J 

The  procedure  used  to  calculate  the  sampled  volume  of  aerosol  of 
type  s  is  as  follows: 

a)  polynomial  representation  of  V  as  a  function  of  A  in  order  to 

determine  V,  and  VA(-.^S  y 

b)  polynomial  representation  of  A  as  a  function  of  V  in  order  to 
perform  the  integration; 

c)  identification  of  cases; 

d)  determination  of  s  ; 

c 

e)  evaluation  of  sampled  volume. 

3.3  Estimation  of  Inertial  Effects 


In  the  previous  sections,  the  interactive  time  of  a  particle 
with  the  air  stream  was  not  considered.  It  was  assumed  that  a  particle 
inside  the  interaction  zone  as  defined  previously  was  automatically 
collected.  But  it  is  now  known  that  a  particle  at  its  critical  posi¬ 
tion  will  never  be  collected.  Hence,  how  far  from  its  critical  posi¬ 
tion  does  a  particle  have  to  be  in  order  to  be  collected  during  the 
sampling  period?  In  order  to  answer  this  question,  the  time  taken  by 
particles  of  s  equal  to  16,  100  and  200  g  um^/cm-  as  a  function  of 
initial  position  R  *  0  and  Z  «  X  BC(s)  with  X’0.1  -  0.99j  ,  and  BC(s), 
the  critical  position  on  the  Z-axis  has  been  calculated  using  [  lj  ,  [  2]  , 
[  3]  and  [ 4 j  . 
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Figure  10  shows  the  results  of  these  calculations.  The  vertical 
lines  represent  critical  positions.  Also  shown  on  the  graph  Is  the 
time  taken  by  the  air  at  a  position  Z  and  R  =  0  to  cross  the  filter 
plane  (Z  =  0).  These  curves  were  obtained  for  a  sampling  flow  rate  of 
50  L/min.  For  a  60-s  sampling  period  it  appears  that  for  s  greater 
than  100,  a  particle  at  0.99  BC  will  be  collected  in  less  than  20  s. 
However,  for  s  =  16,  particles  at  a  position  above  0.94  BC  will  not  be 
collected.  In  order  to  attain  0.99  BC  for  s  =  16,  the  sampling  period 
would  have  to  be  extended  to  150  s.  It  can  be  seen  that  the  Inertial 
effects  are  not  very  important.  However  it  is  important  to  remember 
that  neglecting  them  will  result  in  an  overestimation  of  the  volume 
sampled,  especially  for  small  particles. 


-30  -25  -20  -15  -10  -5  t 

Z-AJ1S  (at) 

FIGURE  10  -  Capture  time  as  a  function  of  distance  along  Z-axis  for 
three  values  of  S 
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4.0  SAMPLING  EFFICIENCY 


The  sampling  efficiency  E  is  defined  as  the  ratio  of  the  sampled 
volume  of  particles  (V.  (  )  +  V  (  ))  of  parameter  value  s  to  the  sampled 

D  3 

volume  of  air  QT : 


E(s) 


V 


)  +  v  ) 

TJT - 


[16] 


Figures  11  and  12  show  the  contribution  to  the  total  efficiency 
(+++)  of  particles  above  ( - )  and  below  ( - )  the  plane  Z  *  0  at  sam¬ 

pling  flow  rates  of  50  and  30  L/min.  Other  parameters  were  h  =  100  cm, 
T  «  60  s  and  rim  size  ■*  2.3  cm.  It  is  clearly  shown  that  for  small  s, 
the  main  contribution  to  the  total  efficiency  comes  from  below  the 
plane  Z  *  0.  As  s  tends  toward  zero,  the  contribution  from  below  the 
plane  Z  =  0  tends  toward  unity;  the  contribution  from  above  the  plane 
tends  toward  zero.  The  contribution  to  E  for  particles  under  the  plane 
is  zero  for  s  with  XC  <  L* 


Figures  13  to  17  show  the  influence  of  height  (h),  sampling 
period  (T),  sampling  flow  rate  (Q)  and  rim  size  (L  -  a). 

The  sampling  flow  rate  effect  on  sampling  efficiency  is  shown 
for  the  two  flow  rates  of  30  and  50  L/min  in  Fig.  13.  Other  parameter 
values  are  T  *»  60  s,  h  °  100  cm  and  rim  size  «  2.3  cm.  As  expected, 
the  sampling  efficiency  increases  with  increasing  sampling  flow  rate. 

Figure  14  shows  the  height  dependence  of  sampling  efficiency  for 

h  «  0,  50  and  100  cm.  Other  parameter  values  were  Q  -  30  L/min, 

rim  size  »  2.3  cm  and  T  -  60  s.  The  difference  in  E  is  quite  large 

between  h  ■  0  cm  and  the  other  two  heights.  The  difference  is  however 

quite  small  between  h  =  50  and  100  cm.  There,  s  was  determined  by 

c 

[l5]  as  66.5  and  133.8  g  pm2/cm3  respectively.  This  is  why  the  differ¬ 
ence  appears  for  values  of  s  greater  than  66.5  g  pm2 /cm3.  For  values 
of  s  greater  than  135  g  pm2 /cm3,  the  volume  collected  comes  only  from 
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®  5C  IK  IK  r  J5f  3K  K  W 
*  (g  u»*/cn5) 

FIGURE  17  -  Influence  of  sampler  rim  size  on  sampling  efficiency 


below  the  filter.  This  is  because  critical  positions  AC(s)  for 
s  >  135  g  um2/cm3  are  below  9  cm,  which  is  the  radius  (L)  of  the  sam¬ 
pler  (6.7  +  2.3  cm) . 

The  sampling  period  (T)  dependence  is  shown  for  three  sampling 

periods  (30,  60  and  120  s)  in  Fig.  15.  Other  parameter  values  were 

Q  =  30  L/min,  h  =  100  cm,  rim  size  =2.3  cm.  The  s^  are  267,  133.8  and 

66.8  g  ura2/cm3  respectively.  There  is  no  significant  difference 

between  these  three  curves.  As  with  the  height  dependence  curves, 

these  curves  cross  one  another  at  s  =  66  g  pm2/cm3.  The  contribution 

c 

of  particles  below  the  filter  becomes  less  important  as  T  increases. 

The  effect  of  rim  size  on  sampling  efficiency  is  shown  in 
Figs.  16  and  17.  Parameter  values  were  h  =  100  cm,  T  •  60  s,  rim 
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sizes  =  0,  2.3,  5.3  cm.  In  Fig.  16,  Q  *>  30  L/mln  and  in  Fig.  17, 

50  L/min.  The  larger  the  rim,  the  lower  the  contribution  from  above 
the  filter.  For  the  limiting  case  of  a  rim  size  larger  than  A  (which 
is  where  the  potential  line  delimiting  the  sampled  volume  crosses  the 
sampler  active  surface  plane),  only  particles  below  the  filter  will 
contribute  to  the  efficiency.  As  the  rim  gets  smaller,  the  contribu¬ 
tion  of  particles  above  the  filter  to  total  efficiency  becomes 
more  and  more  important.  However,  for  particles  with  large  s 
(s  >  200  g  pn^/cm3  ,  Fig.  16;  s  >  340  g  pm2 /cm3  ,  Fig.  17),  only  parti¬ 
cles  directly  below  the  filter  will  be  collected.  This  is  not  due  to  a 
short  sampling  period  but  rather  to  the  relatively  low  speed  of  air  for 
regions  not  directly  in  front  of  the  filter.  Also  shown  in  Figs.  16 
and  17  are  the  curves  of  sampling  efficiency  as  predicted  by  Walton 
(Ref.  2),  which  states  that  for  a  very  thin  rim,  efficiency  is  given  by 
the  expression 

1-/  [17] 

where 

V  :  settling  velocity  of  particles 
6 

V:  airspeed  at  the  filter. 

This  theory  does  not  take  into  account  the  height  of  the  aerosol 
above  the  sampler  and  the  sampling  period,  and  it  assumes  a  very  small 
rim  size. 

5.0  CONCLUSIONS 

A  computer  simulation  has  been  used  to  calculate  the  sampling 
efficiency  of  an  inverted  sampler  with  a  large  entrance  aperture.  It 
has  been  shown  that  the  sampling  efficiency  can  be  represented  by  the 
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sum  of  the  sampling  efficiencies  of  particles  above  and  below  the  sam¬ 
pler.  As  particle  sizes  tend  towards  zero,  the  main  contribution  to 
total  efficiency  comes  from  particles  below  the  sampler.  For  very 
small  particles,  the  sampling  efficiency  tends  toward  unity. 

The  height  and  the  sampling  period  of  aerosols  above  the  sampler 
have  a  strong  Influence  on  sampling  efficiency.  The  sampling  effi¬ 
ciency  of  large  particles  increases  with  height  and  sampling  period. 

However,  the  sampling  efficiency  is  height-independent  for  heights 
2gPr2aT 

greater  than  — - - ,  where  r  is  the  cut-off  particle  size  collected 

9p  a 

for  particles  above  the  sampler  plane.  Note  the  linear  relationship 
between  the  height  and  the  sampling  period  T.  As  expected,  a  higher 
sampling  flow  rate  allows  collection  of  larger  particles.  Finally,  a 
smaller  rim  size  will  provide  a  higher  sampling  efficiency,  especially 
for  large  particles. 

The  approach  is  general  enough  to  be  applied  to  other  inverted 
samplers.  An  improvement  to  this  simulation  would  be  the  use  of  the 
exact  flow-field  calculations  around  the  sampler. 
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APPENDIX  A 


Potential  and  Flow  Lines  Created  by  a  Flat  Disk 
with  Uniform  Velocity  across  its  Surface _ 


The  Laplace  equation  in  cylindrical  coordinates  is: 

32$  S2^  1  3$_  1_  S2^ 

3  22  3  R2  R  3  R  R2  302 

In  the  case  of  symmetry  about  the  Z-axis,  there  is  no 

)c  z 

G-dependence.  A  particular  solution  is  then  4>  =  e  JQ(kR)  for  z  <  0 
where  k  is  a  positive  Integration  constant  in  the  Four ier- Bessel 
space . 


Using  Four ier-Bessel  expansion  and  the  boundary  condition  that 

3  $ 

on  the  surface  (z  =  0),  -  *  f(R),  a  general  solution  is  obtained 

(Ref.  Al),  I .  e . : 


4  =  J  ekZ  JQ(kR)dR  J  f(R’)  J  (kR’)  R’  dR’  [  A- 2 ] 
o  o 

If  the  speed  is  uniform  over  the  disk  of  radius  R  =  a  at  z  =  0, 
and  since 


/  J q( kR ' )  R’  dR'  -  |  Jj  (ka)  [ A- 3 ] 

o 

the  solution  becomes 

QO 

4  =  aV  I  ekz  J  (kR)  J  j (ka)  [  A-4] 

The  expression  for  the  current  lines  is  obtained  using  the 
theorem  of  nul  divergence  for  a  harmonic  function,  i.e.: 


1  3  (pV_)  ,  3  V 
Tr  r  +  7z  Z 


R 


[A-5] 


V  • 


0 
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„  1  3</  ...  13* 

VR  “  "  TT  ST  and  Vz  =  R  3"r 


[  A“6' 


and  furthermore: 


v*  =  e 

v®  3T  r  Jz  z 


v  $  +  v  e 

R  R  2  Z 


[A-7] 


Then,  using  [  A- 5  and  [  A— 6 ] ,  we  obtain 


13*  3*  #  1  3*  3* 

~  IT  5T  '  ST  ’  R  51 X  *  TZ 


Integrating , 


rod)  (  3  $ 

<J<  “  -R  /  5-5-  <*z;  or  *  =  /  _  RdR 


using  equations  A-7  and  A-4  ,  [  A— 4 ]  becomes 


«  kZ 


*  =  -aV  /  fL—JjCkR)  Jj(ka)  dK 


'  A~9] 


Using  the  expression  relating  the  sampling  flow  rate  Q  to  the 

area  of  the  disk,  V  ■=  —9.  ,  the  product  aV  can  be  replaced  by  — — . 

it  a2  1,3 
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APPFNTIX  B 


Note  on  Velocitv  Calculation 


Equation  A-8  has  been  used  to  calculate  current  lines  in  the 

space  R-Z.  A  matrix  of  m  x  n  elements  represents  this  region  of  space 

(m  =  n  =  21).  Velocities  V  and  V  are  calculated  using  the  following 

R  Z 

expressions: 


V  =  -  I 
R  R  TZ 


VRN  (I*J)  = 


*(i+i,J)  - 
Sz 


Z  R  Tr 


*  vzs  (I'J) 


♦  O.J+1)  -  *  (I,J) 
Ir 


where  AZ  and  AR  are  distances  between  nodes  on  the  R  and  Z  axes; 

V  (I,. I)  and  V  (I,J)  are  velocities  on  nodes  along  the  R  and  Z  axes. 
RN  ZN 

AZ  and  AR  have  been  set  to  1.5  and  1.68  cm,  which  is  a/4. 


The  position  of  nodes  for  the  velocities  matrix  is  summarized  in 
Table  B-I. 


TABLE  B-I 


Node  on  R-axis 


Node  on  Z-axis 


(J-l)AR 


i  ( 21-1 )AZ 


J(2J-1)aR 


(I-l)AZ 
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The  velocity  matrix  has  (m  -  1)  x  (n  -  1)  elements  and  they  are 
not  superimposed  one  on  another.  Figure  B1  shows  the  position  of  nodes 
for  ip  with  (*),  with  (X)  and  V  with  (+). 

Internodal  velocities  are  calculated  using  interpolation 
equations  (see  Fig.  B2). 

VR  -  <1-«VR>  <i+1'j» 

+  +  4  W1+1'J+1» 

I  B-3] 

vz  -  (1-4Rvz>  <0-«vz>  VZN(I,J)  +  tzvz  VZN(I+1,J)) 

+  4Rv2  <<1-4Zvz)  V2(I,J*1>  h-  SZvz  Vzk(I+1,J+1)) 

[B-4] 

where  internodal  distances  AR  and  AZ  have  been  normalized  to  1. 
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